Membrane morphology induced by anisotropic proteins 
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There are a great many proteins that localize to and collectively generate curvature in biological 
fluid membranes. We study changes in the topology of fluid membranes due to the presence of highly 
anisotropic, curvature-inducing proteins. Generically, we find a surprisingly rich phase diagram 
with phases of both positive and negative Gaussian curvature. As a concrete example modeled 
on experiments, we find that a lamellar phase in a negative Gaussian curvature regime exhibits a 
propensity to form screw dislocations of definite burgers scalar but of both chirality. The induced 
curvature depends strongly on the membrane rigidity, suggesting membrane composition can be a 
factor regulating membrane sculpting to to curvature-inducing proteins. 



To form the structures of many cellular organelles from 
the Golgi apparatus to the endoplasmic reticulum to the 
mitochondrial inner membrane, lipid bilayers must be 
molded and shaped by an army of proteins jl[ . Among 
these are proteins known to localize to and, collectively, 
induce curvature in fluid membranes [l|-|3[. A flurry of 
theory [^-Q , simulations [^, Q , and experiments [Ifl, |ll[ 
has established that proteins can, in principal, induce 
large-scale shape changes in membranes. Nevertheless, 
hovif these proteins control morphology and what other 
factors are important in determining morphology is still 
poorly understood. 

In this paper, we develop a mean-field model to de- 
scribe how anisotropic, curvature-inducing proteins in- 
duce topological changes in membranes. Using our gen- 
eral model, we show that the magnitude and sign of the 
induced membrane curvature depends on the protein in- 
trinsic curvatures and concentrations as well as on the 
rigidities of the underlying membrane. Since the effective 
membrane rigidity can depend on lipid composition, our 
model qualitatively explains the observed role of mem- 
brane composition in determining curvature 10|, illj . For 
concreteness, we explicitly apply our model to the forma- 
tion of bicontinuous structures in a lamellar phase doped 
with saddle-forming proteins. Note only have bicontin- 
uous, cubic phases been observed in model membranes 
using a variety of proteins [ll|, ll2| , but they are also seen 



in vitro in the inner membranes of the mitochondria of 
starved ameoba [l3j. 

To each bound protein, we associate a position r on 
a two-dimensional surface by projecting the protein cen- 
ter along the normal to the midsurface. A unit vector 
u, tangent to the surface, points along its backbone (see 
Fig. [T]). The second fundamental form of the midsurface 
can be expressed in terms of two principal curvatures ci 



and C2 as hi 



ciBi iSi j -\- c2e2.iG2.j1 where the e^ are 



unit vectors along the principal curvature directions. Fi- 
nally, we define a unit vector t = N x u, where N is 
the surface normal, that points transverse to the protein 
backbone. It will prove convenient to define the angle 
cos 6 = u • ei , where the dot product is taken with re- 
spect to the surface metric, and dimensionless prescribed 



curvatures ^^ = Ci/\cL\- Assuming Hookean elasticity for 
the protein-membrane interaction, we obtain the energy 



E±{v,6 



= ^(CiCos2 + C2sin2 0Tl)' 



+Xx[{^i-^2)cos9smeT^x] (1) 



where the sign refers to the sign of £,l- The first term 
describes the interaction along u while the last term de- 
scribes the transverse interaction, which induces the cur- 
vature S.TCL = ct along t. One can think of the Xi 
as the energetic cost, in units of UbT /2, of binding an 
intrinsically-curved protein to a flat membrane. When 
At = AjT = Ax do we recover the previously studied 
models [I,!]. 

Rather than fixing the protein density, we calculate in 
a fixed chemical potential ensemble. Since the proteins 
on each leaf need not be in equilibrium, we must intro- 
duce different chemical potentials, /i±. If a protein, when 
bound to one leaf of a bilayer, induces a curvature S,i, it 
must induce a curvature —S,i when bound to the oppo- 
site leaf (Fig. [T]). Therefore, we associate /i± with the 
positive or negative sign of Eq. ([Ij respectively. Since 
the proteins only interact through the membrane curva- 
ture, the partition function for either membrane leaf is 




FIG. 1: Schematic model of proteins binding to a bilayer 
membrane. The proteins are at the same position, r, un- 
der a normal projection to the midsurface (dashed) with u 
and t tangent to the membrane. If proteins bound to the 
top layer induce a local curvature +cl along the midsurface, 
an identical protein bound to the bottom surface induces a 
curvature —cl- 



found simply by exponentiating the single protein parti- 
tion function, 



n± 



exp 



dAde 



Ao 



-z±e 



-E±{e,T)/{kBT) 



(2) 



where dA is the area measure on the membrane 
midsurface, Aq is a characteristic cross-sectional area 
of the protein and z± = exp[mu±/{kBT)] are 
the protein fugacities. The average protein con- 
centration and orientation is, therefore, (p(r, 0)) = 
z± exp{-S±(r, 0)/{kBT)}/{2TrAo). 

Finally, we combine Eq. ([2]) with the Helfrich energy 
for a membrane, finding 



knT 



= / dA 



2 (^ - ^o)' 



^if-E^5± 



(3) 



where g±{r) = J d9 {cxp[-E±{9,r)/{kBT)] - 
expl-Ep/ikBT)]}, H{r) = (ci -f C2)/2 is the mean 
curvature, Hq the spontaneous curvature, K{r) = ciC2 
the Gaussian curvature. We have also subtracted off 
a contribution in .g±(r), unphysical in a fixed area 
ensemble, involving the interaction energy between the 
proteins and a flat membrane. It is straightforward to 
generalize Eq. (jS]) to include multiple species of proteins. 

We seek to minimize Eq. (jSj with respect to the mem- 
brane shape. Before proceeding, it is worth noting that 
a straightforward minimization will not capture the cor- 
related fluctuations of the membrane and, by proxy, the 
correlated fluctuations of the proteins themselves. Hence, 
the proteins interact only through the large-scale equilib- 
rium deformations of the membrane. This restriction can 
be corrected by computing the fluctuation corrections to 
Eq. dni). 

In order to establish our main points with minimal 
algebraic complication, we will mostly specialize to the 
extremely anisotropic case Xx = Xt = 0- As in the 
lamellar experiments, we will assume that proteins can 
bind to either leaf with equal fugacity, so that z = 
z+ — z^. Later, we will reintroduce the remaining cou- 
plings to see how this basic picture is modifled. With 
this simplification, Ep = e"^^/^ and E±{r,e)/{kBT) = 
Xl (Ci cos2 9 + ^2 sin^ 61 T l) /2 in Eq. 0. 

As a first step toward a topological phase diagram, we 
ask what combinations of membrane curvatures d and 
^2 minimize the free energy density. Fig. [2] shows a 
typical phase diagram as a function of Xl and fugacity 
z. We generically find three morphologies: a flat phase 
(Ci = ^2 = 0), a saddle phase (^ = — ^2 ~ 1), and 
a spherical phase (^ = ^2 ~ !)• When S = 0, the 
transition from flat to minimal surface occurs along the 
line Xl — 1; therefore, we use R — —0.05k in Fig. [2] to 
highlight the role of R. 

We can understand the transition from the flat to 
curved phases in terms of the two "topological mod- 
uli", K+ = k/2 + R and k_ = —R [ij, |l5[. When 




kAoc| 



FIG. 2: (color online) Longitudinal coupling strength, \l, vs. 
normalized fugacity, z/{kAoCi^), phase diagram for a single 
species with equal binding to both leaves and Ax = At = 0. 
We've chosen k — —0.05k. The dashed line is the stability 
line of Eq. ((5| between the fiat phase (red/dark gray) and 
the minimal surface phase (blue/gray). The spherical phase 
is shown in yellow, with a transition from the K < phase 
that is first order. 



K^- < 0, the membrane becomes unstable to topolog- 
ical rearrangements that induce positive Gaussian cur- 
vature, for example to spherical vesicles [15|. Similarly, 
when K- < 0, the membrane becomes unstable toward 
the formation of negative Gaussian curvature. By ex- 
panding Eq. ([3]) in powers of the curvature, which is 
valid only in the flat phase, we find effective moduli are 
shifted by Ak = 7r(z+ + z-)/iAocl)e-^^/^3{l - Xl)Xl 
and Ak = 7r(z+ + z-)/ {2Aoc\)e-^^/'^XL {Xl - 1) 1] so 
that K+ < when 
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Eq. dS]) predicts the boundary between the flat and 
saddle phases in Fig. [5] quite well. If k < — k/3, 
the transition from the fiat phase proceeds directly 
to spheres, as can also be seen from Eqs. ^ and 
([5]). In contrast, an instability with respect to contin- 
uous perturbations occurs at a critical protein fugac- 
ity (z+ -I- z^)/{AoKcl) > e^^^'^/lSnXLiXL - 1)] above 
which the effective bending modulus k + Ak < 0. There 



K = 




KAgci 



FIG. 3: (color online) Transverse coupling strength, At, vs. 
renormalized fugacity z/(kAocJ^) phase diagram (k = 0) for 
a single species of protein with Ax = Ai, = 10 and with 
equal binding to both bilayer leaves. The solid, black lines 
indicate a discontinuous transition while the dashed, black 
line is continuous. There is a change of symmetry across each 
of these black lines. Cylinders {K = 0) lie on the red line. 
The color intensity is indicates the magnitude of the Gaussian 
curvature. 



is also a shift in the spontaneous curvature of AHq — 
2ti{z+ — Z-) / {AQCL)e~^^ /'^ Xl / {n + Ak) which vanishes 
when the fugacities z+ = z_ = z are balanced. 

Why do anisotropic couplings lead to such rich phase 
diagrams? When \x = At = 0, the protein curvature 
can be accommodated by both K > Q and K < Q since 
they need only find a single direction in which the mem- 
brane curvature is commensurate with the protein's in- 
trinsic curvature. However, the orientational entropy of 
the protein is maximized at umbilics (^i — ^2) while 
the bending energy is minimized at minimal surfaces 
(^1 = —^2)- The dependence on the combination z/k 
arises directly from this competition. 

Terms higher than second order in the curvature aris- 
ing in Eq. ^ are responsible for stabilizing the highly 
curved phases in Fig. [2 These terms typically com- 
pete with potential higher order terms in the Helfrich 
energy. These Helfrich terms are t ypi cally ignored since, 
in the absence of other physics [l6|, they are associ- 
ated with microscopic length scales, i. The higher or- 
der terms in our model, on the other hand, scale with 
{z+ + z_)e~^^/^/(^oc|), the average protein density on 
a flat membrane. If we denote R as the average separa- 



tion between proteins on the membrane, we obtain that 
the regime of validity of our model, with respect to higher 
order corrections, R -C l/{tc\) for our model to be valid 
without these additional corrections. 

If we introduce the coupling Xx 7^ 0, the spheri- 
cal phase increases its range of stability. Even when 
Xx — Xl, however, all three phases remain at higher 
couplings, though at small couplings (Ax — Xl ~ 3) 
the spherical phase preempts the saddle phase entirely. 
When At ^ as well, additional phases can appear. 
If ^T — 0, for example, a "cylindrical" phase, with 
l^il 7^ IC2I and 1^2! ^ 1, is introduced. As an exam- 
ple, a cross-section of the phase diagram along the slice 
Xl = Ax = 10 is plotted in Fig. [H The transitions meet 
at a single critical point at which ^i = ^2 = 0. For com- 
pletely isotropic couplings (At = Ax = Xl), there is a 
"cylindrical" phase (with a small but nonzero K) above 
a critical fugacity and coupling strength. At small but 
isotropic couplings, however, protein entropy dominates 
so that only the flat and spherical phase remain. 

Additional complications arise when we make the tran- 
sition from local phase diagrams, such as those in Figs. 
[2] and [3l to equilibrium membrane morphologies. We 
focus on the case that Ax = At = and k = 0. In 
the case of the spherical phase, the lowest energy state 
is one of monodisperse spherical vesicles of radius c^ . 
In the saddle phase, however, the lowest energy state 
occurs when ^1 = —^2 ~ 1- This is impossible in a 
three dimensional Euclidean space - the membrane en- 
ergy is frustrated by the constraints of geometry. As a 
prototype of a topological transition driven by proteins, 
we consider the lamellar to bicontinuous phase transi- 
tion observed experimentally in the presence of certain 
curvature-inducing peptides 



ll| . Generically, lamellar 



phases that have a tendancy toward negative Gaussian 
curvature can display a variety of complex morphologies 



17| . in many cases, however, these complex layered struc- 



tures can be built from a superposition of simple defects 
[l8|. In fact, screw dislocations are the building blocks 
of more complex minimal surfaces |19j . 

Here, we will consider the free energy of inserting a 
single screw dislocations into the lamellar order. Once 
it becomes energetically favorable to insert dislocations, 
we should expect them to proliferate in the ground state. 
Since screw dislocations of either sign will be degener- 
ate, we should expect the formation of a globally achiral 
structure such as that explored in Ref. pj|]. In terms of 
coordinates {x, y) on a flat reference layer, a screw dislo- 
cation can be described by the multivalued height func- 
tion of a hehcoid, h{x,y) = btan'~^{y/x)/{2n), where 
b = nd is the Burgers scalar, n an integer, and d the 
layer spacing. In the absence of proteins, this is an ex- 
tremum of both the bending and compression energies 
[20|; since we are only after the stability of the lamel- 
lar phase, we will assume the same height function with 
proteins present. Screw dislocations have a large elastic 



6 = 3 




with proteins from a spherical to saddle phase. Should 
curvature-inducing proteins take advantage of our mech- 
anism for sculpting membranes, tuning the membrane 
moduli by adding cosurfacants or changing lipid mix- 
tures will induce a corresponding transition in membrane 
Gaussian curvature. 

In summary, we have studied topological transitions in 
membranes induced by anisotropic, curvature-inducing 
proteins. We predict a surprising variety of morpholo- 
gies having both positive and negative Gaussian curva- 
tures depending on the protein densities and membrane 
moduli. In a regime of negative Gaussian curvature, this 
induces a transition from a lamellar phase to one with 
screw dislocations of both chiralities in the ground state. 

This work was funded by the National Science Foun- 
dation through DMR-0846582. 



FIG. 4: (color online) Longitudinal coupling strength, \l, 
vs. renormalized fugacity z£_ /{Bd Aocl) critical line above 
which helicoids are energetically favorable for R = Q. The flat 
phase is stable in the region A_l < 1. Color (shading) is added 
as a guide to the eye. 



contribution to their line tension, arising from the devi- 
ation from the equilibrium layer spacing near the core, 
of the form r = Bb'^ / {2^11'^^'^), where B is the bulk 
modulus and ^ is a microscopic core size [20|. Due to 
the divergence as ^ -> 0, we neglect an additional core 
energy in comparison to the elastic line tension. 

Since 7J = 0, this line tension must be balanced with 
the contribution from Eq. ([3]) , which we evaluate numer- 
ically. We find that a screw dislocation becomes energet- 
ically favorable when compared to the flat phase above 
the critical line shown in Fig. |4l Increasing fugacity 
prefers an increasing burgers scalar, with both chiralities 
degenerate. This critical line occurs at larger Al than 
predicted by the local phase diagram. The transition 
also depends strongly on the bulk modulus of the lamellar 
phase, though we note that the layer spacing is typically 
set by d ~ ^J{K/d)/ B 20] . If we use ^ ~ d for the core 
size we obtain that z^^/(i?(i'^AoCi) ~ zcLd/{K,AQc\), so 
the horizontal axis of Fig. |4] is roughly comparable to 
that of our local phase diagrams when cl ^ 1/d. 

Finally, the dependence of our phase diagrams on 
membrane rigidity imply that the induced morphology 
also depends on lipid composition, since lipids with in- 
trinsic curvature can modify the membrane rigidity [Tf. 
For an anisotropic lipid with A^ < 1, for example, our 
results indicate that Ak > 0, suggesting that these 
intrinsically-curved lipids can drive a membrane coated 
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